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In this paper we propose and extensively stndy mimetic f{G) modified gravity models, with various 
scenarios of cosmological evolution, with or without extra matter fluids. The easiest formulation is 
based on the use of Lagrange multiplier constraint. In certain versions of this theory, it is possible 
to realize accelerated expansion of the Universe or even unified evolution which includes inflation 
with dark energy, and at the same time in the same theoretical framework, dark matter is described 
by the theory. This is achieved by the re-parametrization of the metric tensor, which introduces a 
new degree of freedom in the cosmological equations and leads to appearance of the mimetic “dark 
matter” component. In the context of mimetic f(G) theory, we also provide some quite general 
reconstruction schemes, which enable us to find which f{G) gravity generates a specific cosmological 
evolution. In addition, we also provide the general reconstruction technique for Lagrange multiplier 
f{G) gravity. All our results are accompanied by illustrative examples, with special emphasis on 
bouncing cosmologies. 


I. INTRODUCTION 

The recent observational data [Il,i , have increased the scientific interest in cosmology, since early-time acceleration 
is favored by present observations. In addition to this early-time acceleration, in the late 90’s, the late-time acceleration 
(including that from modified gravity Si) was also verified by the observations S. Since then a great task for 
cosmologists is to consistently describe these two acceleration eras, and if possible, to include them in the same 
theoretical framework as it was first proposed in Q. Modified gravity S, is one appealing candidate for a complete 
cosmological theory, in the context of which late-time acceleration can be consistently described but more importantly, 
the unification of early-time acceleration with late-time acceleration S 0 can be achieved (for a general review, see 
™). There are many types of modified theories of gravity with higher-derivative terms, such as f{R) modified gravity 
f{G) theories of gravity [l^, etc. In this paper, we shall introduce a new type of mimetic f{G) theory of gravity, 
which as we shall demonstrate, can successfully describe a variety of cosmological scenarios. It is inspired by mimetic 
dark matter models [lol - [T^ and serves as a generalization of both f{G) and mimetic dark matter models. 

The mimetic dark matter models [lol - fl^ . are generalizations of general relativity, at which general relativity is 
modified in a minimal way, so that the conformal symmetry not to be violated. In the theoretical formalism of 
mimetic dark matter theories, an extra internal degree of freedom appears, which can describe dark matter via such 
mimetic gravity theory 0 (see also Refs. [Hllil)- In order for the conformal symmetry not to be violated, the metric 
is not regarded as the fundamental variable of this gravitational theory, but is decomposed into a new metric and a 
scalar field, with the latter being the extra gravitational degree of freedom which actually describes dark matter. 

A recent study in which the formalism of mimetic gravity was used in the context of f{R) gravity was done in Ref. 
[l3l |. As was shown in that work, in the context of mimetic f{R) gravity and specifically of modified f{R) gravity 
with Lagrange multiplier constraint, it is possible to realize the unification of early-time and late-time acceleration 
eras as well as having Dark Matter in unified scenario. 

It is of great importance to understand how universal is such an approach to the unification of the three dark sectors 
of the Universe. 

In this paper we consider modified Gauss-Bonnet f{G) gravity, which was introduced as gravitational alternative 
to dark energy in Ref. [l^ . but generalized in such a way that the Lagrangian multiplier constraint is used, as done 
in Ref. [l5l| . The specific choice of the Lagrange multiplier constraint, renders such a theory as a mimetic modified 
Gauss-Bonnet gravity of certain type. The accelerating cosmology generated by such a mimetic gravity, is studied in 
detail and we also investi gate which f{G) gravity can realize various cosmological scenarios, with special emphasis to 
bouncing cosmologies [l6l - l^ . We need to note that the idea of obtaining a bouncing solution in the frame of Lagrange 
multiplier modified gravity was studied for the first time in [^ . which is very relevant in spirit to our approach. In 
addition, we also introduce the mimetic gravity formalism, in the context of a general Jordan frame f{G) gravity. As 
a general remark, we must note that within the context of mimetic f{G) gravity and Lagrange multiplier modified 
f{G) gravity, it is possible that there exists a large number of f{G) theories which can consistently describe a given 
cosmological evolution. 
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This paper is organized as follows: In section II we present the theoretical formalism of mimetic Gauss-Bonnet 
gravity. The most important argument is based on the fact that the metric is redefined, and all cosmological equations 
are obtained by varying the redefined action with respect to redefined metric = —g^'^dp(j)dn4>9iiv and mimetic scalar 
field (j) instead of In section III we investigate how various cosmological scenarios can be realized in the context 
of mimetic f{G) gravity, with special emphasis to cosmological bounces. In section IV, we consider a generalization of 
f{G) theory, based on the introduction of a Lagrange multiplier, which includes the mimetic case. Several examples 
are worked out and also comparison to the standard Einstein-Hilbert gravity with Lagrange multiplier, worked out 
in [llj l is performed. Moreover, we also compare the resulting modified f{G) gravity, mimetic or Lagrange multiplier 
modified, with the standard f{G) gravity picture. The conclusions follow in the end of the paper. 


II. MIMETIC f{g) GRAVITY 


Consider a general f{G) gravity, with the Jordan frame action being equal to. 


5 ' = 






( 1 ) 


where Sm stands for the action of matter-fluids, which induces the energy momentum tensor T^i, in the field equations 
(for some incomplete list of works on cosmology of such theory, see M)- We adopt the standard geometrical 
conventions of Einstein-Hilbert gravity, with regards to the commutative connection of the underlying Riemannian 
spacetime, R is the Ricci scalar, and the function f{G) corresponds to a generic globally differentiable function of the 
Gauss-Bonnet topological invariant G, which is equal to, 

G = R^- + Rp^xaR^''^^ , (2) 


where R^^, and RpuXa are the Ricci and Riemann tensors, respectively. We adopted the signature for the Riemannian 
metric as (— h -I--I-) and also we assume that = SttG/c'*’ = I, where G is the standard Newtonian gravitational 
coupling. 

In order to obtain the mimetic f{G), the following parametrization of metric is assumed [TMl: 

= -g^'^dp(t>da(l)gpu- (3) 

By varying the metric we obtain the following relation. 


Sg^l.l^ = g^^Sgru,g‘^''dp(j)da(l)gp„ - g^'^dp(t>da(l)dgpy- 


-2gP'^dp(l>da5(t)gpy. 

In addition, by varying the action o, with respect to the redefined metric g^y, instead of the standard Jordan frame 
metric g^v, and with respect to the mimetic scalar cj), the resulting field equations take the following form. 


R 





( 4 ) 


-b8 


R 1 

R^pua T Rpvgap Rpagvp. Rp.v>gap “t” Rpagvp H ^ igpi^go-p g^ag^p^ fg -t- {fgG f{G)) gpy~\~ 


-\-dp(j)d,j(j) ( —R J- 8 ( —Rpcr + —Rgpa 


VPV^fg+4{fgG-f{G))] = 


= Tpy + dp(t)dv(j)T, 

where fg stands for fg = df{G)/dG- We also note that the covariant derivative V^, acts on vectors as V^I4 = 
dpVi, — ^p^Vx and accordingly it acts on the metric g^,^. 

By varying the action o with respect to the scalar field (j), we obtain the following equation, 

(^-Rp, + ^Rgp,^ V^V^fg + 4{fgG - fiG)^ - T)^ = 0. 


-R + 8 


( 5 ) 
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Assuming a spatially-flat Friedmann-Robertson-Walker (FRW) metric, with metric 

= —dt^ + a^{t){dx^ + dy"^ + dz^) 

the scalar curvature and the Gauss-Bonnet invariant have the following form: 


R = 6 





( 6 ) 


g = - 24 ^ = -2AH^(H 4 - H"^). 

where H{t) denotes as usual the Hubble rate H{t) = a{t)/a{t). From Eq. ([3]) we find that 

gf"’'df,4>d^(j) = -1 (7) 

and since the scalar field (j) depends only on the cosmic time, this results to the constraint tp = t. Taking into account 
the aforementioned constraint, we can obtain the (t, t) component of the expression given in (|3|): 

2H + 3R2 + iqh{H + _ f(^g)) = _p. (8) 

The same equation results if we consider the (r, r) component. By integrating Eq. ®, we obtain the following 
relation. 


- R + 8 (^-Rp^ + ^Rgp?j + ^{hG - fiG)) + p-3p = -^ (9) 

and finally this results to, 

H + 2H^ + + 4R (2H + 3H^) ^ + ^{fgg - /(S)) + f " | (10) 

where C is an arbitrary constant. 

Combining Eqs. ([5]) and m, we obtain: 

H + + 4H(2H - = -i (, + p) - (11) 

where for simplicity we redefined the constant C. It is convenient to introduce the function g{t), which is defined to 
be equal to, 


which satisfies the equation: 


9{t) = 


dt ’ 


+ 4Hi2H - H^)g{t) =-H-\{p + p) - (12) 

at 2 

Eor (7 = 0 Eqs. ([5|) and (fTTl) coincide with the Friedmann equations in ordinary f{g) gravity. If the right hand side of 
Eq. (fT3] is zero, that is, B{t) = —H —\{p + p) — ^ = 0 then the solution to this equation can be written as follows, 

= 5o ■ (13) 

In the above equation, g^ is a constant and Hq = H{0). In the case that the right hand of Eq. (1121) is non-zero, that 
is when B{t) ^ 0, the solution gii) reads, 

(^) Jo H{T)dT^ . 


9{t) = 90 


(14) 
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Therefore one can specify the cosmic evolution of the Universe, in terms of the Hubble rate H{t) and then obtain the 
function g{t). Then by integrating this, with respect to the cosmic time, we can obtain the function /(t), that is, 

/eW = J 9{t)dt. 

Using the explicit dependence of the Gauss-Bonnet invariant on the cosmic time t, one can get the explicit form of 
the function t[Q) and therefore find the function fg{Q). Then, upon integrating with respect to Q, we easily obtain 
the mimetic f{G) gravity. In the next section, we shall present some illustrative examples, in order to demonstrate 
explicitly how the reconstruction method works, and also in order to see the differences between Jordan frame mimetic 
f{G) gravity and ordinary Jordan frame f{G) gravity. 


III. COSMOLOGICAL MODELS WITH VARIOUS EVOLUTIONS FROM MIMETIC f(g) 

A. A general example 


We start off with the cosmological evolution of the Universe, with a Hubble rate equal to. 

Ho 


H{t) = 


1 + at’ 


a > 0 


which could correspond for example to early-time acceleration from the moment t = 0 to some later time t, with 
t >> a~^. The corresponding scale factor for the Hubble parameter (IHI A() is equal to, 


a{t) = oo (1 J- 


(15) 


The corresponding Gauss-Bonnet invariant is, 


G = 2A 


H^{a-Ho) 


(1 + atY 

The solution of Eq. (II2I) with zero right hand side {B = 0), is equal to, 

git) = 50 (1 + . 

For an empty universe without matter and mimetic dark matter {p = p = C = 0) the function g{t) is equal to, 

\Ho/a +2 1 1 + at 


git) = 5o (1 + at) 

Then, the resulting expression for the function /q(t) is equal to. 


4iJo a + Hq 


8aHo a + Hq 


Using the Gauss-Bonnet invariant (IHI All , we can express the function fgit), as a function of the Gauss-Bonnet 
invariant G, 

fg = jgg-V3^ 


D = - 


SaHoia + Hq) 


{2AH^{a - Ho)) 


2/3 


Finally, integrating the above equation, with respect to the Gauss-Bonnet invariant G, we obtain the resulting non- 
mimetic fiG) gravity. 




( 16 ) 
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For a Universe filled with matter and mimetic dark matter, which means that C ^ 0, p = poa p = 0, the 
corresponding solution g{t) has the following form, 


g{t) = 50 (1 + 


1 1 + at 

4:Ho a + Hq 


G 1 
dFTo - a ml 


(1 + , 


C = C + po, 


(17) 


and the corresponding solution fg{t) reads, 


fait) = 


(18) 


where E stands for. 


E = 


C 


a 


ml (4i7o - a)(4a - Ho) 


{2ml{a - Ho)) 


4/3-Ho/a 


The resulting mimetic /(!/) gravity easily follows, 


f{g) = Fog-^°/^ + 3 Dg^/^ + ^g-l/3+go/a_ (jg) 

SHo — a 

By looking Eqs. m and (HSl), we can see the difference between non-mimetic and mimetic f{g) gravity. Notice that 
in the mimetic gravity case, a given cosmological evolution can be realized by a family of f{g) functions. For example 
by putting C = —po and therefore the last term can be omitted in Eq. (ED. It is worthy to note that in ordinary 
f{g) gravity with matter fluids, we have no such freedom in choosing the f{g) function, and a unique solution for 
f{g) exists (in (fT71) one need to change C —>■ po)- 

Let us study another case, in which the cosmological evolution is described by the following Hubble rate. 


H = 


Hq 

1 -k /32t2 


We consider the simple case for which 50 = 0 and p = p = C = Q, and the solution g{t) is equal to. 


9{t) 


+ pH^)il + Hot - /3H^) 
2Ho{Hl + 4 / 32 ) 


The Gauss-Bonnet invariant is. 


2AHl{2pH - Ho) 
+ 1)3 


For t —^ 00 we have the following asymptotic behavior: 


( 20 ) 


... 

~ 2Ho{Hl + 


g~ 


A&Hl 

PH^ ■ 


Therefore, the function fg{g) is equal to. 


fg{G) 


^ 14/5 


2{48Ho)-^/mo{Hl + 4py 
and upon integrating with respect to g, we finally obtain the /(!/) gravity, 

5^14/5 


,- 4/5 


fiG) 


2{A8Ho)-^/^Ho{Hl 


■m 


^ 1 / 5 , 


( 21 ) 
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B. Reconstruction of mimetic f{G) models from bounce cosmological evolution 

Using the method we described in the previous sections, it is possible to reconstruct an f{G) model for given 
cosmological evolution with scale factor a(t). In this section we shall consider bouncing cosmologies, and we investigate 
which f{G) model can generate such a cosmology, in the presence of matter and mimetic dark matter. 


1. Superbounce from mimetic f{G) gravity 

In the context of f{G) modified gravity, various bouncing Universe scenarios can be materialized straightforwardly, 
which in the context of Einstein-Hilbert general relativity could only be realized if the null energ y co ndition is violated 
for a period of time. For reviews and an important stream of papers on bouncing cosmology see [l6l-l20j and references 
therein. A bounce evolution has the following characteristics: The Universe starts with a contraction up to a point 
where a non-zero minimal radius is reached, where the scale factor is non-vanishing too, and then the Universe starts 
to expand again. In this way the singularity problem of the scale factor, which is known as the initial singularity, is 
resolved in a self-consistent way. Therefore, bouncing cosmology serves as a convenient way to avoid singularities of 
crushing type, where the geodesics incompleteness occurs. Let the bounce occur at some time U, then for t = the 
Hubble rate is zero, that is H{ts) = 0, while during the expanding phase, H{t) > 0, and during the contracting phase, 
H{t) < 0. 

Therefore the bouncing cosmology is one appealing candidate for the consistent description of the Universe’s evo¬ 
lution. Within this context, we shall make use of two bouncing cosmologies for our reconstruction examples. 

We start of with the superbounce model [2l|, which has the following scale factor. 


a{t) = {-t + ts)' 


( 22 ) 


while the corresponding Hubble rate is, 




c 2 (-f -h U) 

Using Eqs. (I21D and (1^51) and substituting to Eq. (HI, the solution g{t) easily follows, 


(23) 


git) = ^c^goH^it - t,)^i-t + Uf - - Q^i-t + 


-2 


-k 


c^Ci—t -k tsY 
8 — -k 


c^it-ts) ip + p) 


2 (2 -k 3c2 -k c4) 


The Gauss-Bonnet invariant for the Hubble rate of Eq. ((23ll is equal to, 

g = -96c2t^(2c-k4c2t^) 
and by solving Eq. (l25l) with respect to t we obtain two real solutions, 


ti — 


-2y23i/^c4yG (c® (-2 -k c^) G) ' -k c®Gt 


2 ^23^/^0470 (c® (-2 -k c2) G) ' -k c®GG 

- (26) 


By substituting Eq. (l26l) to Eq. (l24l) and integrating over G, we obtain two different functional forms of mimetic fiG) 
gravity, which are of the following form, 


/(a) = -12(-2 + c2) itiG)-ts) 
2c^Ci-t{G)+tsf-^ 


-1+^ / 4 (Sgpgg -k (1 - 4goiJg)) 
° V (_2 + c2)" (2 + 3c2) 

c*{tiG) - G)^(P + P) 


(27) 


(-4 -k c2)4 (8 - 14c 2 -k 5c4) (2 -k c^) (2 -k 7c2 -k bc^) 

where t{G) = 0,2 with 0,^2 being defined in Eq. (E51) . 

Let us here compare the Jordan frame mimetic fiG) modified gravity solutions we found, that generate the su¬ 
perbounce (1221) . with the solutions of non-mimetic fiG) gravity in the Jordan frame. The non-mimetic fiG) gravity 
solutions that generate the superbounce (l2^ are the following M- 


tlG - tlG-2^i-n + a)a^Gil + G) 

j\y) = - 77—772 -> fiy) = — 


4a -k 4a2 


4a(l -k a)G 


( 28 ) 
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where a = 2/c^. Comparing Eqs. (I27p and (l28l) . we can see the obvious functional differences. But there is also 
another important difference, having to do with the freedom in choosing the solutions in the case of mimetic f{G) 
gravity, an issue we already discussed earlier. Notice also that there is a similarity between the mimetic and non- 
mimetic f{G) gravities, which is that in both cases, the resulting f{G) gravity consists of two distinct solutions, when 
the superbounce cosmological evolution is taken into account. 


2. Generic bounce scenario from mimetic f{G) gravity 


Another interesting bouncing cosmology scenario is described by the following scale factor, 

a{t) = 


(29) 


with c some arbitrary constant and the scale factor is normalized again so that at the bouncing point, which is t = 0 , 
the scale factor is equal to a(0) = 1. The corresponding Hubble rate is particularly simple, and it is equal to. 


H{t) = 2ct 

Using the Hubble rate (l30l) and by substituting to Eq. Cl, we get the function g{t), 

Ce-2ct^ ect^p ^ct^p C^eTi[V2^ct\ v^erfi(v^t) 


sW = 


64c^t 16c^t 


4c^t 


64c^t 64c^t 


32c5/2 


16 c 3/2 


(30) 

(31) 


go77oV^erfi(v^t) i/7rperfi (i/ct) 


4c3/2 
'1. 3 . 
^ 2 ’ 2 ’ 

16c^ 


+ 


64c5/2 
rl. 3 . 
<2J_2 ’ 

8 c 


ct 2 F 2 {^-, ct'^) t 2 ^ 2 (^; f;ct^) pt 2 F 2 {F ct'^) ^^ 2 ^ 2 ( 5 ; i;ct^) 

32^^ 


64c3/2 

A- 3. 

^ 2 ’ 2 ’ 

32c2 


where the functions erfi(a;) and erf(x) denote the imaginary error function and the error function respectively, while 
the function pFq(a\ 6 ; z) stands for the generalized hypergeometric function. The Gauss-Bonnet invariant in this case 
is equal to. 


g = -96c^t^ (2c + 4ch^) 

which can be solved with respect to t, and we have the following positive solution. 


t=\ — 


4c 


V24c6 - c^G 
8^c4 


(32) 


(33) 


Having this at hand, and by substituting to Eq. (ED), we easily obtain the function fg{G)- Then by integrating with 
respect to G-, we can get the resulting fg{G) function, but we omit the resulting expression, since it’s functional form 
is too large and complicated. Note however that all the integrals can be evaluated analytically, as the reader may 
convince himself. 

Before closing, we study here a final example, for which we assume that go = 0, in which case the function g(t) is 
equal to. 


git) = - 


a\t) 

4d2(t) 




(34) 


Consider a cosmolo gica l evolution with scale factor a{t) = (1 -b at)'*'. The case a < 0, 7 < 0 corresponds to the 
Big Rip singularity [2^, which occurs at tg = 1/a. The accelerated expansion corresponds to a > 0, 7 > 1. The 
corresponding solution for g{t) reads. 


1 


5 W = ((1 + - (1 + <^tf) - 


C 


8720 ^( 47 —!) 

The corresponding Gauss-Bonnet invariant G is, 

e = 24^)“'‘3 -J’’ 

(1 -b at)4 

so the resulting expression for the f{G) gravity reads, 

f{g) = Ag^!^ -b -b 


^(1-b at)^+^ - (1-b , 7 7^1/4. 


(35) 


(36) 



























IV. f{g) GRAVITY IN FORMULATION WITH LAGRANGE MULTIPLIER 


In this section we consider the formulation of /{G) modified gravity by using a Lagrange multiplier, and the impact 
of this multiplier to the cosmological evolution of various models. One can impose a condition on the scalar field by 
adding to the standard action for f{G) gravity the following term, 

+ £/((/.)), 

where A is the Lagrange multiplier and U{cj)) is a function of the scalar field, in principle arbitrarily chosen. We also 
include in the Jordan frame action of the f{G) gravity, the action of a scalar field (canonical or phantom) (j) with 
arbitrarily chosen scalar potential which is of the form: 

50 = y d'^xy/^ -V{4>)). 

The cases e = ±1 describe canonical and phantom scalar field respectively. For e = 0 the conformal symmetry of the 
system is respected and the corresponding gravity with Lagrange multiplier is equivalent to the mimetic f{G) gravity 
considered above (for U = 1 and U = 0). (The case with non-zero V maybe called as extended mimetic theory). The 
full Jordan frame action of the f{G) gravity with Lagrange multiplier and in the presence of the scalar field and other 
matter fluids, is the following, 

S = J (^1 + fiG) - eg>^''d^<t>d^4, - V{cf) + X{g>^^d^cl>d,cj) + !/(<(.)) + (37) 

By varying this action with respect to the metric g^i,, we obtain the following two equations of motion: 

+ 24573^^ + f(g) - fg[g)g = p + 4‘^ + U(^) - A(02 + (38) 

-2H- _ iqh{H + + fgiG)G - fiG) =p + e^^- U(<(.) - A (<^2 _ ( 39 ) 

where we assumed a flat FRW metric. Also, by varying action (EZl), with respect to the Lagrange multiplier, provides 


us with the following constraint for the scalar field: 

<^2 _ ; 7 (^) = 0 . (40) 

Finally for the scalar field we have the equation of motion: 

2dtiiX - e)<j>) + 6H{X - e)(p - V'icj)) + XU'icf) = 0. (41) 

For a Universe dominated by collisionless matter (p = 0), the scalar potential U(t) and the Lagrange multiplier 
function A(t), as function of the cosmic time read, 

Vit) =2H + + £02 + + imiH + - fgiG)G + fiG), (42) 

dt^ dt 

A(t) = fj)-^(P-+H + AHi2H - g 2 ) ^/g(g) 4g2 ^ + e- (43) 

\ 2 dt dt^ ) 


Equation (HD) follows easily by combining Eqs. (1551) and (1551) . Therefore, for a given fiG) one can specify the evolution 
of the scale factor (or the corresponding Hubble parameter) and obtain the scalar field potential which is responsible 
for this evolution. In addition, the inverse process is easy to be realized, that is, for a given scalar potential and for an 
arbitrarily given cosmological evolution, it is possible to reconstruct the corresponding fiG) gravity. In the following 
sections, we exemplify the two aforementioned reconstruction methods. 
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A. Some Examples for f{G) = gravity and (p = y/TJot 

The simplest choice for the function f{G) is the following, 

f{G)=Ag^. (44) 

In this section, we shall assume that the scalar field (p is related to the cosmic time t as, (p = y/Uot, which means 
that 17((/)) = Uq = const. For the bounce case with scale factor given in Eq. (1^ and Hubble rate given in Eq. (1501) . 
the corresponding potential of the scalar field, provided by Eq. (l42l) . takes the following form, 

3 

V{t) = Y. CuiaUo^p^r + Eq, (45) 

n—0 


where the parameter Vq, which is equal to. 


Vo = 3(2c)2c/o-+ 4c + eUo 

stands for the potential of the scalar field in the case of General Relativity (which can be achieved if A = 0), and C„ 
are numerical coefficients. For completeness, we provide the full expression for the scalar potential in the Appendix. 
The corresponding Lagrange multiplier function A(t), which is given in Eq. (I43L is equal to, 


,,,, 36864Ac^ 2c p 221184AcV 98304Ac^(j>^ 

\{(p) = -77-he-'' 


Uo 






(46) 


In the same way, we may easily find the scalar potential and the auxiliary function X{4>), for the case of the superbounce 
with scale factor and Hubble rate given in Eqs. (l22l) and (|23)) respectively. The corresponding scalar potential is. 


VW = 


12288A (-12 + 44c2 - 75c^ + 28c^) 


+ 


c^® {-tsVUo + (p)^ 
c^'^y/Uo {—tsVUo + (pY + c^^(p‘^ — 2.^/Uo (—6 + 2c^ + cHsep)) 


(47) 


C^® {-tsVUo + (pY 

while the Lagrange multiplier function X{(p) equals to, 

24576A (-4 - 120^ + 7c'^) 

X[(p) — 


+ 


(-Gi/C/o + (p) 

c^^ {-ts^/Uo + (pf' {- 2 cHs^/Uo(p (p + epA) + A 4 A [p + ecp"^) + Uo (-2 + Atl [p + ecp'^))) 


(48) 


+ (p) 

Another example that produces Little Rip cosmology, is described by the following Hubble rate, 

H{t) = i7oexp(at), a > 0 . 

with the corresponding scalar potential being of the following form, 

3 

v{(p) = AHq exp( 57 (()) Y Cna^~'^HQ eyipiri'ycp) + Vo, 

n—Q 


(49) 


Vo = 3 Hq eyip{2j(p) + 2Hoa exp^'jcp) + eU, 7 = aU 

In order to better understand what is the difference between the Lagrange multiplier mimetic matter with a scalar 
potential, which was studied in Ref. [TTI |. and Lagrange multiplier f{g) gravity, we shall study in detail some examples 
that were presented in [HI- It is worth recalling some features of Lagrange multiplier mimetic matter, with scalar 
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potential, and for a detailed presentation, the reader is referred to 0. The action of the mimetic matter with scalar 
potential and with Lagrange multiplier, is given below [^ . 


S = 


dx^ 


--i? + A - 1) - 1A((^) + 


(50) 


with Lm denoting as usual the matter Lagrangian. In the context of the theory developed in [m, the Lagrange 
multiplier is equal to. 


A=i(G-r-4I/(<^)), 


(51) 


with G and T standing for the trace of the Einstein tensor and for the trace of the energy momentum tensor 
respectively. The form of the Lagrange multiplier appearing in Eq. (1511) is obviously different from the Lagrange 
multiplier appearing in Eq. (l4^ . but one can easily verify that if we take t/(())) = 1 and also f{Q) = 0, these become 
identical. We now investigate the new features that the f{G) modification brings along. We choose again the f{G) 
gravity to be of the form given in Eq. (|4^ and also we assume that (p = y/Uot. We shall start our comparison analysis, 
by studying the inflaton case presented in [HI, where it was shown that a cosmological evolution with scale factor 
equal to. 


a{t) = (52) 

with a a constant parameter. The potential that generates this kind of evolution, always in the context of Ref. [HI, 
is equal to. 


V{P) 




(53) 


In contrast, in the case of the Lagrange multiplier f{G) modified gravity, the potential that generates the cosmological 
evolution described by Eq. dSlD, is equal to, 


V{P) 


2^ f— 002 256 ^ 0 * 5/202 lO88y4a304 14O8^o5'/206 64^0^0® 

-^ + VGoe+ —+ + 9731/3 


(54) 


which is clearly much different than the one in Eq. (15311 . In addition, the corresponding Lagrange multiplier function 
A(0) is equal to. 


^ p 128^o®/2 (9V3C/2 - 2717ov/S02 - 2v^o04) 

^ +- 2 TO 5 - <“> 

which is obviously different from the corresponding one in ordinary Lagrange multiplier mimetic matter theory given 
in Eq. (1511) . It is obvious that within the context of Lagrange multiplier f{G) modified theory of gravity, there are 
many ways of generating various cosmological scenarios, since there is much freedom in choosing the potential t/(0) 
and also the function f{G)- In this way, almost any cosmological scenario can be generated, by suitably choosing the 
two aforementioned functions. Before closing this section, we shall study another example presented in [TTj |. which 
describes mimetic matter as quintessence. The scale factor of the corresponding cosmological evolution, is equal to 
0 . 

a(t) = (56) 

which actually describes the behavior of mimetic matter when the Universe is dominated by some other matter fluid 
present, with equation of state p = wp. The potential that generates the cosmological evolution (l55|) . is equal to [HI, 

U(0) ^ ^ (57) 

where a is some arbitrary constant. In the case of Lagrange multiplier f{G) gravity, the potential that can generate 
the cosmological evolution of Eq. (l56l) . is equal to, 

f4096A(l + 3u;)(67 + u;(149 + 84w)) - 972^(1^+^)°-^° + 729(i+yfe0^ 

U( 0 ) = ^ 0 u. 


729(1 + ^)808 


( 58 ) 
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and the corresponding Lagrange multiplier function X{4>), is equal to, 

8192AU^{l + 3w)i2^ + 2lw) 2Uo p 

~ ^ ^ 729(1 + u;)70io 3(1+ w)^^ ^ Ip' 


(59) 


which are clearly different from the ones given in Eqs. (IS3 and m- 

As we already mentioned, in the context of Lagrange multiplier f{G) modified gravity, there is much freedom in 
reconstructing the scalar potential for the scalar field, since the function U{(j)) and the function f{G) itself can be 
arbitrarily chosen. Let us consider as a final example the case for which the function that relates the scalar field cj) to 
the cosmic time t is of the following form, (j) = •yZTolnt, so the function U{<j)) is of the following form, 


U{(j)) = C/oexp(-2f/p 


(60) 


Assuming that. 


Hq = at 

and also that the function f{G) is again f{G) = AG^ , the corresponding scalar field potential V{(j)), is of the form, 

3 

V (</)) = Aa^ exp(2UQ ^^^(/)) ^ Cna" exp(n[/Q ^^^(()) + Vb, (61) 

n—0 


where Vq in this case stands for, 

Vq = 2a + 3a^ exp(217g'^^^^) + ello exp(—217g'^^^^). 

In the same way, more elaborated universe evolution unifying inflation, dark energy and dark matter maybe presented. 


B. Reconstruction of f{Q) gravity given the Lagrange multiplier function and the scale factor 


In this section we consider the inverse reconstruction method, which enables us to find the Lagrange multiplier 
f{G) modified gravity, given the scale factor of the cosmological evolution and the Lagrange multiplier function. Let 
us explicitly demonstrate how this method works. We start off with Eq. (H51) . which can be rewritten as: 


(A(,)-.W=-f-A-4a(,)d(^f 1^0, 


Solving this equation with respect to fg (t) we have: 


(A(t2) - e)'^^(i2) - H{t2) - 


2a^{t2) J 


(62) 


(63) 


Note that for A = 1, Eq. (|6^ in fact coincides with (|35ll (with (7 = 0) for g{t) = dfg/dt, as was expected. The 
potential of scalar field as a function of time is cast as follows. 


V{t) = (2A(t) - e)(^^(t) + 3H‘^(t) - p{t) + 6a{t)H{t) [ ( (A(ti) - - H{ti) - 


i{ti) 


2a?[ti) 


(64) 


-fg{t)Q{t) + f{i)- 

We assume that the function that relates the scalar field (f> and the cosmic time t is denoted as t We also 

write the scale factor a{t) in the following form. 


a{t) = exp(/i(t)). 


(65) 


Then, the Hubble parameter is simply written in terms of h(t), as H(t) = h(t). Considering for simplicity the case 
that no matter fluids are present, that is p = 0, the potential of the scalar field given in Eq. (l64)) . becomes, 

V{(j)) = 3h^{y{(j))) - - 6e^^'^^‘>‘^h{y{cj))) J dcj)iy'{(j)i)e-'^^y^^^^^ + (66) 
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-fg{(j})G{y{(t))) + J d(l)ifg{(l)i)y'{(t)i)gy{y{(j)i)) 

and therefore, the function fg as a function of the scalar field is equal to, 

fgi^) = -\J + Hym)^ + (67) 

Then, for a given cosmological evolution, in terms of some specific scale factor, by choosing the arbitrary Lagrange 
multiplier function A(())), we can easily obtain the corresponding /{G) gravity, and also finally find the resulting form 
of the potential V{(/)), which we used earlier for the derivation of the f{G) gravity. 

Let us exemplify how this inverse reconstruction method works, by using some illustrative examples. For simplicity, 
we shall assume that the field (f> is related to the cosmic time as, ^ = \/Uot, so we define the function y((/)) to be the 
following, 


y(0) = (68) 

which we extensively use in the following considerations. Of course, one can in principle choose a different function, 
but this will only perplex the equations, without changing the qualitative feature of our reconstruction method. 

We start off with the superbounce, with scale factor and Hubble rate appearing in Eqs. (USD and (USD respectively. 
The function h{t) is in this case equal to. 


h{t) = In {—t + ts)'^ 

In addition, we take the function A(t), to be equal to, 


\(t) = ct 

Then, by combining Eqs. dSTD, dSHl), dsni) and GOl), we obtain the function fg{4>), which is equal to, 

. (l+c^)tab^ (- 3 +c^)ta 0 ^ (- 2 +c^) 0 ^ 


c^l- 


h{(t>) = 




vdTo 




2VUo 


5VUo 


32 (2 - 3c2 + c4) Uo 


^(j){—2ts + (jd) ^2c^ (2 + c^) tie — 2(? (2 + (?) tgecj) + C/q ^—8 + ^ + 2? ^2 + 


Uq 


64 (-4 + ?) Uo 


(69) 


(70) 


(71) 


Then, the f{G) gravity easily follows, by using the expression for the Gauss-Bonnet invariant, which for the case of 
the superbounce is given in Eq. (1^51) . by integrating with respect to G- The final expression is quite large and for 
simplicity reasons we give it in the Appendix. In the same way, by using Eqs. (IMl), dMl), dMl) and dZOl), we obtain 
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the potential of the scalar field, which is, 


n<^) = -7r + 

U(\ 


12 


2c(j) 6 c — 2 (/)^ + (—+ (/)^)) 


12 


c^{ts-4’Y (2-3c2 + c4)C/o/^(t,-0) 

(2 + c^) t^e — 2c^ (2 + c^) tse<p + Uq ^—4 + —h 2c^ ^1 + ^ 


c 2 (-4 + c4)17o(t;, - (t>Y 


(72) 


4 _ 2 \ („7 ( c^tt<P , (l+c")t^0^ _ _ {-3+c^)ts4>^ (-2+c")0" \ 

{ts —0)^ y ^ ^ \/C/o \/t/o ■'/t/o 2-\/Uo 5-\/T7o y 

c4(t^ - <)))232 (2 - 3c2 + c4) [/o 

( c4(t,-,^)2 - e2(t,-0)2 ) (c'‘(/)(-2ts + ()i) (^2 c 2 (2 + c^) tfe - 2 c 2 (2 + c2) + Uq (-8 

c‘^{ts - (f))^64 (-4 + c^) Uo 

3 (-2 + c2) (^-2cHl + hcHl + icHl - 40t7o/^ - ^cHly/U^e - ^c^y/Uo (4[/o + tie) + IQc^y/U^ (6C/o + tfe))) 

5c^C/y^ (4 — 4c2 — + c®) (ts — 

I 40C/o^^ , 8c^<#. , 20c^(cts-\/t77e) ln[-t„+</.] 

(2+c2)(t3-0)2 + -l+c2 ' -2+c2 

5c4C/y' 


96 ( 
H- 



In the same way, we can easily obtain the f{G) gravity for the case of the bounce with scale factor (1291) . Using the 
same conventions as in the case of the superbounce, the resulting expression for the f{G) gravity is. 


4C/2''\ 


4 ^16c®-^ 


(73) 


( v^l 6 c® - ^ 3 ( 4 


3/2 


16'^ c 


y/TTUo(2cUo + e) 


xG 22 


- 1 , 0,1 

34 I 0,0,-2,-i 


12c3Uo - Vl44c® - 6c4GC/o 


48c3 


1 / —36c37r + 3-\/144c® — bc^GTr 

4 V ? 




12c3Go - Vl44c® - 6c4GGo 


48c3 


with G™", the Meijer G-function. The resulting expression for the scalar potential is too complicated, so we provide 
the details for it in the Appendix. 

As a final example, we shall consider a non-monotonic evolution for the Hubble rate which has the following form, 

2 2 
H{t) = - 


t tf-t 

Notice that a Big Rip singularity occurs at t = tf. It easily follows that. 


h(4.) = ln[-J-lny^ C„ 


h{(j)) = 


2yM , SyUo 


= —IT 


2Uo 


Without getting into too much detail, the function fg(4>) is equal to. 


/eW = - 


16Uocj)l 






(t^f = \/Uotf. 


5 ^ ^ 4 ^ 3 


25 


4>] 2 


+ <^2 + 


(74) 
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If,. A-2(y(<^i))eMy(0i)) „-/.to(0.)) 


+17s,»(0.) 

while the corresponding scalar field potential is equal to, 
Uo 


d4>2e 


K4>2) 
y'i(p2)J ' 


n<^) = 


{(j)f - (j)f 


+ I2(t>)(t)-^ + e{ll(t>) + 3(j)}(j) - 15(/>/(/)2 + + 24(;!)^(/)ln 0) - 60Ci(?!>/(/)) + (75) 


2/ (w J \y wi 


Notice that the last integral can be evaluated analytically, but we omit the resulting expression for the sake of brevity. 
In conclusion, with this inverse reconstruction method, if the Lagrange multiplier function A(0) and the scale factor 
a{t) are specified, the corresponding Lagrange multiplier f{G) gravity follows, as we explicitly demonstrated in the 
previous examples. Of course, it is to be understood that if the integrals cannot be evaluated analytically, further 
assumptions should be taken into account, but the general theoretical framework is the same. 


V. DISCUSSION 

It is worth discussing certain issues related to the mimetic f{Q) gravity models we studied in this paper. In particular 
one possible question that comes to mind when dealing with this kind of theories is the following: all the models we 
constructed in the present study involve an f{G) function, with G being the Gauss-Bonnet term which in the usual 
case is harmless, since it can be integrated out topologically. Now the question is if the theory constructed by using 
a function of the Gauss-Bonnet invariant, contains any extra degrees of freedom or not. In our case however, note 
that the f{G) gravity is considered here as an effective string-inspired gravity theory. Due to non-linear functional 
dependence on G, the function f{G) is not a topological invariant anymore. The introduction of the f{G) term, 
introduces to the theory an effective scalar degree of freedom, as in many cases the model maybe rewritten as one 
containing a string-inspired V{(j))G term. In addition, in the context of the mimetic f{G) gravity, there is an extra 
scalar degree of freedom, that of the scalar field 4>, which accounts for the dark matter content of the theory. Now we 
come to another important issue that we need to briefly discuss. As we demonstrated in this paper, the mimetic f{G) 
cosmological scenario can serve as a theoretical framework that at the same time can describe inflation, dark energy 
and even dark matter. Although it is quite appealing to have a geometric description of the cosmic history of our 
Universe, an important question rises, with regards to the inflationary particle content of the theory. In particular, 
for inflation, the most important and successful achievement was to predict a nearly scale invariant power spectrum 
which is verified by CMB experiments in the past decade. However, this achievement is based on the fundamental 
assumption that inflation is realized by a slow-rolling scalar field in the early universe. Without such a scalar field, 
how can the CMB spectrum be realized in modified gravity models? 

This question is very interesting and we shall try to briefly address this, so let us first consider the standard mimetic 
gravity, presented in |1^ . ini , where the authors considered the standard Einstein-Hilbert mimetic extension. As it 
was shown in [TTj |. the metric perturbation is related to the higher derivative of the scalar field’s linear perturbation 
(5(/), and in effect, the resulting picture is that, the perturbations behave as dust, with the speed of sound being zero, 
even for mimetic matter. Therefore, the quantum fluctuations of mimetic matter are ill-defined, and therefore the 
large scale structure cannot be explained by the minimal Einstein-Hilbert mimetic matter theory. An extra scalar 
field is required to play the role of the curvaton or explicitly modify the mimetic dark matter model, in order for 
the large scale structure being consistently described. Therefore, the mimetic f{G) could potentially provide such an 
extension, by offering extra degrees of freedom to the theory. The corresponding extra scalar degree of freedom we 
discussed earlier, may play such a role in the resulting theory. In [llj |. the authors modified by hand the minimally 
extended mimetic Einstein-Hilbert Lagrangian, thus being able to produce a red-tilted spectral index. It would be 
therefore quite interesting to investigate whether this extra scalar term can originate by an f{G) theory or from any 
other mimetic extension of modified gravity. However, this task exceeds the purposes of this paper and we hope to 
address such issues in the future. 
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VI. CONCLUSIONS 

In this paper we considered the generalization of ordinary Jordan frame f{G) gravity, to mimetic f{G) gravity. 
In this new theory, in principle it is possible to resolve two important cosmological problems, firstly to obtain the 
Universe’s accelerated expansion without the need to introduce extra helds (inflaton, dark energy scalar) and secondly 
to provide an answer to the problem of dark matter. Specifically, within the context of mimetic f{G) gravity, the 
cosmological equations of motion are almost the same with these that result in ordinary Jordan frame f{G) gravity, 
with the only difference being that the contribution of mimetic matter appears in the mimetic f{G) gravity case. 
After having provided the formalism of mimetic f{G) gravity, we were able to realize various cosmological evolutions 
with mimetic f{G) gravity, and highlighted the difference between ordinary and mimetic f{G) gravity. In principle, in 
the case of mimetic f{G) gravity, there is much more freedom in providing the f{G) gravity that can realize a certain 
cosmological evolution. This can be one of the attributes of the theory, since a quite large number of cosmological 
evolutions can be successfully realized in the context of mimetic f{G) gravity. We also compared the resulting mimetic 
f{G) gravities, to standard mimetic gravity models that exist in the literature. 

We also modified the ordinary Jordan frame f{G) gravity, to include a Lagrange multiplier term. After describ¬ 
ing in detail how the cosmological equations are modified within this new formalism, we provided two quite general 
reconstruction techniques, which in principle can be useful for the realization of various cosmological models. Par¬ 
ticularly, it is possible in the context of Lagrange multiplier f{G) gravity, given the cosmological evolution and the 
f{G) gravity, to find the scalar field potential V{(f>) and the Lagrange multiplier function X{t). In addition, the inverse 
procedure is possible, that is, given the cosmological evolution and the Lagrange multiplier function A(t), it is possible 
to reconstruct the corresponding f{G) gravity. We supported our theoretical considerations by using some illustrative 
examples, for which we applied both reconstruction methods. 

As a general remark, we must note that in the context of both mimetic f{G) gravity and Lagrange multiplier f{G) 
gravity, there exist a large class of f{G) gravities that can realize a specific cosmological evolution. We believe that 
this feature is an attribute of the theory, since it is possible, by suitably choosing the parameters and functions, to 
get analytic results. This was not always possible in ordinary Jordan frame f{G) gravity. 

An issue we did not address in this paper is the realization of singular cosmology, in the context of mimetic f{G) 

S vity or Lagrange multiplier f{G) gravity. By singular cosmology we mean the appearance of finite time singularities 
(with the Big Rip being the most elaborated one) in the cosmological evolution. Among all the finite time 
singularities, the most interesting ones are the non-crushing types singularities, for example. Type II or the mildest 
among them, the Type IV [2^ . . Since these singularities can be consistently incorporated in scalar-tensor theories 

[ 2 ^ . it would be interesting to realize such singular cosmological evolution from mimetic f{G) and from Lagrange 
multiplier f{G)^ and compare the results to ordinary f{G) gravity. Special emphasis for this task, should be given 
near the Type IV singularity of course. Work is in progress towards this direction and we hope to report on this issue 
soon. 

Finally, let us discuss another issue having to do with the possible observational distinction between the various 
modified gravity theories. Indeed, these theories can in principle be very degenerate with regards to their observational 
prediction so the fundamental question that comes to mind is, whether there is any possible way of distinguishing 
these theories at a pragmatic level. One good idea is to these these modified gravity models by using weak leasing 
data, owing to the fact that the trajectories of light may differ, depending on different geometric terms. For some 
relevant works on this issue see for example [ 2 ^. In addition, the growth index may provide some useful insights on 
that respect, see for example [30l| . 


Acknowledgments 

This work is supported in part by project I4-02-3II00 (RFBR, Russia) and project 2058 (MES, Russia) (AVA), by 
MINECO (Spain), projects FIS20I0-I5640 and FIS20I3-4488I (SDO). 

Appendix: Resnlting expressions for f{Q) gravities and for scalar potentials 

Here we quote the complicated expressions we mentioned in the main text of this paper. The full expression for 
the potential appearing in Eq. (1451) is equal to. 


V{(l)) = 4c -f - 73728AcH^ - 626688Ac®t^ - 540672Ac^t® - I47456Ac®t® -f y/u^e 


(76) 
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In addition, by integrating Eq. (ED, with respect to t/, we obtain the f{G) function, which is equal to, 
^ + 20(4-4c2-c4 + c 6 )(-t, + ^) 


(77) 


X tl(2cHl^/ui^(lO - 20^/u^ + 15Uo - + 5c^tl(^- 4 + 6^/u^ 

- AUo + C/q/^) - 4o( - + cHl ( - 10 + 10\/Uo - + 2U'tj 

- - 4 + 6v^- 4C/o + e - 5 c4( - 2+ C/q) (2tle - 2tl^Qe + [/q(4 + ) 

+ lOc^^ — 2\/Uq + Uq^ — 2 <g\/Z/oe + C/o ^6 + 

24 ( - 1 + yZio) (2c5t3 ( - 1 + v^) " + chi ( 


3(4-4c=-c<+c«)(-f, + ffl) 

- 4chl ( - 1 + - 8t/o^/" - 2cHI ( - 1 + v^) ^ yhToe 

- 2 c 4 ^/l^ (tie - 2tl v^e + C/q (2 + - 2tl y^e + C/q (s + t 

^ ^ 1 “f \/t7o^ \/Uo + 2c^tg ^ — 1 + \/Uo 


- 1 


C/o) ( 


1 +VC^o 


( 4 - 44 »-c 4 + cs)(-(. + ^^ 

X ( - 1 + 4\/cio) + 4Ul - chl^/Woi^l - 3Uo + 2Ul'‘^') + 1 + y^) Voe 

- Qc^Uq (Uo + tle- 2tl \/4oe + t^t/oe) + c'‘t/o (st^e - yhj^e + t/o (2 + 37 
4c^7s ( “ 1 + a/Ed) C^o (c(^4 - 2tsy/Upj + c^tsy/Uo + v^e - c^a/Edc) 


(2-34» + o4)(-i, + ^) 



Cts 

'3 - 4E!7o) 

+ c®4 

' - 1 + 2y!7o) 

+ VUqc — y/Uoe 

In 

4 , HG)' 






2 - 3c2 + 


J 


where 4>{Q) stands for, 


m = 


or equivalently, 


^{.Q) = 


2y/2?hl^c^yfG (c® (-2 + c2) G) + c®G4 
c®G (EiTo)-' ’ 

-2y23i/4c4EG (c® (-2 + c2) + c®G4 


(78) 


c®G (EtTo) 


-1 


(79) 


In conclusion, we have two different types of f{4>iS)) gravity, which can be found by substituting the two different 
forms of (j){G) to Eq. (iTTl) . Finally, the corresponding potential for the f{G) gravity of Eq. (|73l) . which corresponds 
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to the bounce (l29ll . is equal to, 


^_^ _ 4c(/) i 2 c ^^2 _ 6yce=‘^"v^(2cl7o + e)^erf [^c^] 


Uo 


+ 9Qc^(j)^ (2c + Ac^cj)^) 


Uo 


( 


V^{2cUo + e)G^J 0 o’ -4 ~ 


0,1 


c0 ( 4 (3 + Ac(j)‘^) + 12c^/^^^/Uo{2cUo + e)4’^Gll 


64c2Go 

- 1 , 0,1 
0 , 0 ,- 2,-4 


— C(j)‘ 


3/2 


+ 


C(j) ( Sy/TTVU^{2cUo + e)(/)G23 i g q _ i 


Uo 

0,1 


21 / '-’i 

2 


— C(t>‘ 


u; 


oil 


where G!J1" is the Meijer G-function and erf(a;) is the error function. 


(80) 
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